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In the frame of flavor 5(7(2) Nambu-Jona-Lasinio model with Ua(1) breaking term we found 
that, the structure of two chiral phase transition lines does not exist at low isospin density in real 
world, and the critical isospin chemical potential for pion superfluidity is exactly the pion mass in 
the vacuum. 
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It is generally believed that there exists a rich phase structure of Quantum Chromodynamics (QCD) at finite 



' temperature and baryon density, for instance, the deconfinement process from hadron gas to quark-gluon plasma, the 
transition from chiral symmetry breaking phase to the symmetry restoration phase [1], and the color superconductivity 
[2] at low temperature and high baryon density. Recently, the study on the QCD phase structure is extended 
to finite isospin density. The physical motivation to study isospin spontaneous breaking and the corresponding pion 
' superfluidity is related to the investigation of compact stars, isospin asymmetric nuclear matter and heavy ion collisions 
t-H , at intermediate energies. 

While the perturbation theory of QCD can describe well the properties of the new phases at high temperatures 
and/or high densities, the study on the phase structure at moderate (baryon or isospin) density depends on lattice 
' QCD calculation and effective models with QCD symmetries. While there is not yet precise lattice result at finite 
(SJ • baryon density due to the Fermion sign problem [3], it is in principle no problem to do lattice simulation at finite 
isospin density [4]. It is found [5] that the critical isospin chemical potential for pion condensation is about the pion 
mass in the vacuum, /ij ~ m^. The QCD phase structure at finite isospin density is also investigated in many low 
energy effective models, such as chiral perturbation theory [4,6,7], Nambu-Jona-Lasinio (NJL) model [11-13], ladder 
f--*) , QCD [8], strong coupling lattice QCD [10] and random matrix method [9]. 

One of the models that enables us to see directly how the dynamic mechanisms of chiral symmetry breaking and 
restoration operate is the NJL model [14] applied to quarks [15]. Within this model, one can obtain the hadronic mass 
spectrum and the static properties of mesons remarkably well [15,16], and the chiral phase transition line [15-17] in 
Qh. the temperature and baryon chemical potential (T — ^b) plane is very close to the one calculated with lattice QCD. 
Recently, this model is also used to investigate the color superconductivity at moderate baryon density [18-22]. In 
the study at finite isospin density, it is predicted [11] in this model that the chiral phase transition line in T — \lb 
plane splits into two branches. This phenomena is also found in random matrix method [9] and ladder QCD [8]. Since 
the NJL Lagrangian used in [11] does not contain the determinant term which breaks the !7a(1) symmetry and leads 
to reasonable meson mass splitting, it is pointed out [13] that the presence of the L breaking term will cancel 
the structure of the two chiral phase transition lines, if the coupling constant describing the 1/^(1) breaking term 
is large enough. However, this coupling constant is considered as a free parameter and not yet determined in [13]. 
Another problem in the NJL calculation at finite isosipn density is that the critical isospion chemical potential for 
pion condensation /ij = m w is not recovered in the model [12]. In this letter we will focus on these two problems in 
the frame of NJL model with Ua(1) breaking term. We hope to derive the critical isospin chemical potential exactly 
and try to fix the chiral structure by fitting the meson masses in the vacuum. 

The letter is organized as follows. We present at finite temperature and baryon and isospin densities the chiral and 
pion condensates in mean field approximation and meson masses in Random Phase Approximation (RPA) in Section 
2, determine the coupling constant of the Ua(X) breaking term in Section 3, and then analytically prove the relation 
H c j = m n in Section 4. We conclude in Section 5. 

II. NJL MODEL AT FINITE TEMPERATURE AND BARYON AND ISOSPIN DENSITIES 

We start with the flavor SU(2) NJL model defined by 

£ = -tpii^d/j, - m + M7o)^ + £>int , (1) 
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where too is the current quark mass, /j, the chemical potential matrix in flavor space, 

with hb and \xi being the baryon and isospin chemical potential, respectively, and the interaction part includes [13] 
the normal four Fermion couplings corresponding to scalar mesons a, a 0} a + and a_ and pseudoscalar mesons 7r , ir + 
and 7T_ excitations, and the 't-Hooft [23] determinant term for breaking, 
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£int = £ [(<KV>) 2 + OMsTaV') 2 ] + y [det ^(1 + 7 5 )V> + det $(1 - ts)V] 

a=0 

= 1(G + K) [{Wf + (^ 75 tV) 2 ] + \{G - K) [(^) 2 + (^ 75 V0 2 ] • (3) 

For if = and m — 0, the Lagrangian is invariant under C/s(l) ^Za(I) ® Sf/y^) ® SUa(2) transformations, but 
for K ^ 0, the symmetry is reduced to Ub(1) ® SUv(2) ® SUa(2) and the Ua(1) breaking leads to a and a mass 
splitting and ir and 77/ mass splitting. If G = K , we come back to the standard NJL model [15] with only er, 7ro,7r+ 
and 7r_ mesons. 

We introduce the quark condensates 

a u = (uu) , a d = (dd) , (4) 

or equivalently the a and a condensates 

a = (tptp) = (uu + dd) = a u + a d , 

a = (i>T3tp) = (uu - dd) = a u - a d , (5) 



and the pion condensate 



7T - - 1 - 

— = {^i^r+i)) = (i>i^ 5 r-ip) = -^(tpir^tp) , (6) 



where we have chosen the pion condensate to be real. The quark condensate and pion condensate are, respectively, the 
order parameter of chiral phase transition and pion superfluidity. By separating each Lorentz scalar in the Lagrangian 
(3) into the classical condensate and the quantum fluctuation, and keeping only the linear terms in the fluctuations, 
one obtains the Lagrangian in mean field approximation, 



G + K 2 



A„/ = i>S-)i, - G{ai + aj) - 2Ka u a d - , (7) 

where is the inverse of the mean field quark propagator, in momentum space it reads 

0-1 /m_ f 7^ + ^70 ~M U i(G + K)nj 5 \ . . 

m ' lj ~^ t(G + K)ir l5 r^+^o-Md) W 

with the effective quark masses 

M u = mo - 2Ga u - 2K<r d , 

M d = too - 2Ga d - 2Ka u . (9) 

The thermodynamic potential of the system in mean field approximation can be expressed in terms of the effective 
quark propagator, 

G -\- K T 

n(T,VB,»r,(Tu,(Td,n) = G{al + a d :) + 2Ka u a d + ^— tt 2 - - lndet5"}(fc) . (10) 

The condensates a u , o d and 7r as functions of temperature and baryon and isospin chemical potentials are determined 
by the minimum thermodynamic potential, 
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aa u da d dn 



(11) 



It is easy to see from the chemical potential matrix and the quark propagator matrix that for [1b = or fij = the 
gap equations for a u and ad are symmetric, and one has 



o"u = &d = — , a = , 

M u = M d = M q =m a -(G + K)a . 

The quark propagator in flavor space can be formally expressed as 



S m f(k) = 



S U u{k) S ud (k) 
Sdu(k) S d d(k) 



(12) 



(13) 



From the comparison with the definition of the condensates (4) and (6), one can express the condensates in terms of 
the matrix elements of S m f, 



*u = -J 
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7T = 



d 4 fc 
d 4 k 



(27T) 

d 4 k 



TV [iS uu {k)\ , 
4 Tr [iS d d(k)] , 



Tr[{S ud (k)+S du (k))j 5 ] 



(14) 



where the trace is in color and spin space. 

In the self-consistent mean field approximation, it is well known that the meson masses Mm as the bound states of 
the colliding quark- antiquark pairs are determined as the poles of the meson propagators in RPA at zero momentum, 



1 - (G + K)U M (k = M M , k = 0) = 

for a and n, and 

1 - (G - K)U M (ho = Mm, k = 0) = 
for r)l and a, where ITm is the meson polarization function 

^-^Tr [T* M iS mf {jp + k)T M iS mf {jp)] 

and Tm and T* M are the interaction vertexes 
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(18) 



with t± = (n ± iT2)l\pl. Doing the trace in flavor and color spaces, one has 
—J^tr [S uu (p + k)S dd (p)} , 
—^tv[S dd {p + k)S uu {p)] , 



(19) 



n 7r+ (fc) 



-2N r 



d 4 p 



tr [^ 5 S uu (p + k)j 5 S dd (p)} , 
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n 7r _ (fc) = -2N c i J 7^4 tr [l5Sdd(p + k)j 5 S uu (p)] , 

/ d 4 
^-^tr [S uu {p + k)S uu {p) + S ud {p + k)S du {p) + S du (p + k)S ud {p) + S dd (p + k)S dd (p)} , 

IV (fc) =n wo (fc) 

f d 4 p 

= ~N c i J ^2^4 tr bs«5«u(p + k)j 5 S uu (p) - j 5 S ud {p + k)~f 5 S du {p) - -f 5 S du (p + k)j 5 S ud {p) + ~f 5 S dd {p + k)j 5 S dd (p)} , 
now the trace is taken only in spin space. 

III. CHIRAL PHASE STRUCTURE 

We now consider the QCD phase structure below the minimum isospin chemical potential \i\ for pion superfluidity. 
Since the pion condensate is zero, there is only chiral phase structure in this region. 

The simple diagonal matrix <S~^(fc) in this region makes it easy to calculate the matrix elements of the effective 
quark propagator, they can be expressed explicitly as 

s (k)- A + 70 - A " 70 



k -Ex k + E 2 
S dd {k) = t-^-et + 



k -E 3 ko + Ei 1 
S ud (k) = S du (k) = , (20) 

with quasiparticlc energies 

Ei = E u — [i u , E 2 = E u + fi u , E 3 = E d — [i d , E 4 — E d + fj, d , 



E u = Vk 2 + Ml , E d = ^k 2 + Ml , (21) 

and energy projectors 

A M = l/ 1± 7o(7-k + M M )\ (22) 



2 V E Utd 

After performing the Matsubara frequency summation in (14), the gap equations determining the two quark conden- 
sates as functions of temperature and baryon and isospin chemical potentials read 



(2tt) 3 ^fWTMl 
d 3 k M, 

with the Fcrmi-Dirac distribution function 



(2?r)3 7=^ (1 - f{Es) f(E 4 )) , (23) 



/(*) = -^TT ■ ( 24 ) 

We now first discuss the real case, N c — 3. It is well known that in the absence of isospin degree of freedom 
the temperature and baryon density effects on chiral symmetry restoration are different [15]: The chiral condensate 
drops down continuously with increasing temperature, which means a second order phase transition, but jumps down 
suddenly at a critical baryon density, which indicates a first order phase transition. At finite baryon density and finite 
isospin density, while the density behavior of the u- and d-quark condensates are different to each other, they may jump 
down at the same critical baryon chemical potential or at two different critical points. In the case without considering 
the C/a(1) breaking term, the two critical points do exist, and therefore, there are two chiral phase transition lines 
[11] in the temperature and baryon chemical potential plane at fixed isospin chemical potential. What is the effect of 
the C/a(1) breaking term on the QCD phase structure? The calculation in the NJL model showed [13] that if there 
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exists a structure of two chiral phase transition lines depends on the coupling constant K of the determinant term. 
Let 



_ 1 / G-K\ _ K 

a ~2 v'gTk) ~gTk 



(25) 



it was found [13] that at m = 60 MeV the two-line structure disappears for a > 0.11. This means that the two-line 
structure is true only at small ratio a. However, a is a free parameter in [13] and the two coupling constants G and 
K are not yet separately determined. 

Before the discussion of the problem in real world with nonzero current quark mass, it is instructive to analyze the 
phase structure in chiral limit with mo = even though pion condensation will happen at finite m- One can regard 
this as the limit of m — > 0. From the gap equations (23) one can clearly see that when one of the quark condensates 
becomes zero, the other one is forced to be zero for any coupling constants G and K ^ 0. Therefore, there is only 
one chiral phase transition line at any K ^ in chiral limit. The two-line structure happens only at K — 0. In this 
case, the two gap equations decouple, the two chiral phase transition lines are determined by 



1 - 12G 



1 - 12G 



d 3 k 1 
f d 3 k 1 
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/(w-f -f)-/(w + f + f 







['-/(M-(f-«)-f)-/(M + (f-«) + f) 

the difference in baryon chemical potential between the two critical points at fixed T and m is 

Afj, B {T,fj,i) = 3 m . 



= 



(26) 



(27) 



In real world there are four parameters in the NJL model, the current quark mass mo, the three-momentum cutoff 
A, and the two coupling constants G and K. Among them mo, A and the combination G + K can be determined by 
fitting the chiral condensate <r, the pion mass m^ and the pion decay constant f n in the vacuum. In [13] mo = 6 MeV, 



A = 590 MeV, and (G + K)A 2 /2 = 2.435, for which a = 2(-241.5 MeV) 3 



m^ 



140.2 MeV, and f n = 92.6 MeV. To 



determine the two coupling constants separately or the ratio a, one needs to know the r\i- or a-meson properties in 
the vacuum. 

In the vacuum the two quark condensates degenerate, and the pole equations determining the meson masses are 
much simplified as 



4A^ C (G + 



1 - AN C (G - K) 



AN C {G + K) J 



d 3 k 1 E\ 


-Ml 


(27r) 3 E q E 2 - 


-ml/4 


d 3 k 1 E\ 


- M 2 


(27r) 3 E q E% - 


- ^a/ 4 


d 3 k 1 


El 


{2v) 3 E q El- 


-ml/ '4 


d 3 k 1 




WE q E%- 


-m 2 ,/4 



= 



= , 



(28) 



with quark energy E q = wk 2 + M 2 It is easy to see that if K = 0, the masses of o— and a- mesons degenerate and 



the masses of 77'- and 7r-mesons degenerate, and if K = G as in the standard NJL Lagrangian, the a- and ryZ-mesons 
disappear in the model. 

The mass equation for 7r-meson can be used to determine the combination G + K, which is already considered in 
[13], and the mass equation for ^/-meson is related to the combination G — K and then to the ratio 



TT 1 



1 



G + K 12 J 



d 3 k Eg 



(29) 



Combining with the above known parameters mo, A and G + K obtained in [13] and choosing m,, = 958 MeV, we 
have a — 0.29 which is much larger than the critical value 0.11 [13] for the two-line structure at m = 60 MeV. In 
fact, for a wide mass region 540 MeV < m^, < 1190 MeV, we have a > 0.11, and there is no two-line structure. 
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To fully answer the question if the two- line structure exists before the pion condensation happens, we should consider 
the limit fii — [fj where the difference in baryon chemical potential between the two critical points is the maximum, 
if the two-line structure happens. From the approximate result of lattice simulation [5] and the exact result of NJL 
analyze in next section, we know that the critical isospin chemical potential for pion superfluidity is equal to the 
pion mass in the vacuum, /ij = m n . From the analysis above, at this value the isospin chemical potential difference 
between the two critical points corresponding, respectively, to a u — and = in chiral limit and with K = is 
A/ib = 3m w . In real world with ttiq ^ and K ^ 0, with increasing coupling constant K or the ratio a the two 
lines approach to each other and finally coincide at about a = 0.21 which is still less than the value 0.29 calculated 
by fitting m^, = 958 MeV. Therefore, the two-line structure disappears if we choose m v , = 958 MeV. In fact, for 720 
MeV < m v , < 1140 MeV, we have still a > 0.21 and the two phase transition lines are cancelled in this wide mass 
region. When the 77/- meson mass is outside this region, the (1) breaking term is not strong enough to cancel the 
two-line structure, but the two lines are already very close to each other. Considering we have used the 77' mass in 
three flavor world, we should use a smaller 77' mass in two flavor world. Using the approximate relation m 2 , oc Nf one 
can estimate m,< w 780MeV which is in the region where a > 0.21. 

At this stage one may ask in which case the two chiral phase transition lines at finite isospin chemical potential 
appear. To answer this question we pay attention to the relation between the 77' mass and N c . It is believed the 77' 
mass decreases with increasing N c [24,25]. At sufficiently large N c , the Ua{^) breaking induced flavor mixing effect 
can be neglected and there will be two chiral phase transition lines. 



IV. CRITICAL ISOSPIN CHEMICAL POTENTIAL FOR PION CONDENSATION 

Since both thermal motion and nonzero baryon number will increase the critical isospin density for pion superfluidity, 
the minimum isospin chemical potential /uj corresponds to the parameters T — and \ib = 0. In this case the pion 
condensate is zero and the two quark condensates degenerate, the pole equations for the meson masses are then 
reduced to 

d 3 k 1 El -Ml 



r r/ 3 k 1 E 2 - M 2 



1 - AN C (G -K) J 



d 3 k 1 E 2 - M 2 



1 _ 4We(G + /o /|Li._J|_ = , 



< G+K »/(^^-; - 



r E iE^-(M, ± ±^,y/4 

where the quark mass M q satisfies the same gap equation as in the vacuum, 

(2^)3^- =°> ( 31 ) 

namely, 

M q (T = 0, jj,B = 0, f-tj < n c j) = M,(0, 0, 0) . (32) 

Therefore, from the comparison with the meson mass equations (16) in the vacuum, we derive the isospin dependence 
of the meson masses for m < fj,j, 

M CT (/i/) = m a , 

M aa {jil) = m a , 



G 



M a± (fu) = m a± =FM/ > 
M^Jpi) = , 
M n± (fj,r) = m n± T M/ > 

M v ,(^i) = m v , . (33) 

We see that the mesons with zero isospin charge keep their vacuum values, the mesons with positive isospin charge 
drop down linearly in the isospin chemical potential, and the mesons with negative isospin charge go up linearly in the 
isospin chemical potential. Since m a > m n , the above isospin dependence will firstly break down at fij = m n . Beyond 
this value the mass of ir + becomes negative and it makes no sense. This gives us a strong hint that the end point of 
the normal phase without pion condensation or the starting point of the pion superfluidity phase is at ^ = . 

To prove this hint we should consider directly the isospin behavior of the order parameter of the pion superfluidity, 
namely the pion condensate. Deriving the matrix elements S u d and Sdu from (8) and then performing the Matsubara 
frequency summation, the gap equation determining self-consistently the pion condensate as a function of isospin 
chemical potential at T = and \ib = reads, 



.. , d 3 k / 1 1 

l-2N r ■ 



•XG + K) J 



(2tt) 3 



,y/(E q -m/2) 2 + (G + K) 2 n 2 \J (E q + Hi/2) 2 + (G + K) 2 it 2 / 



. (34) 



Obviously, the minimum isospin chemical potential /ij where the pion superfluidity starts is controlled by the square 
bracket with ir = 0, 

r H 3 k 1 E 2 

l-ANXG + K) / — T - ^^ = o. (35) 

>} (2tt) 3 E q E 2 - {^f /4 

From the comparison of (35) with the pion mass equation (16) in the vacuum, we find explicitly that the critical 
isospin chemical potential \x\ for the pion superfluidity phase transition at T = = is exactly the vacuum pion 
mass Wtt, 

H c j = m v . (36) 

We emphasize that this conclusion is independent of the model parameters and the regularization scheme, it is a 
general conclusion in mean field approximation for quarks and RPA approximation for mesons. Unlike the chiral 
phase structure which strongly depends on Ua(^) breaking coupling K and color degree of freedom N c , the critical 
isospin chemical potential pX{ does not depend on them. 



V. CONCLUSIONS 



We have investigated the two flavor NJL model with Ua(X) breaking term at finite isospin density, as well as at 
finite temperature and baryon density. With a self-consistent treatment of quarks in mean field approximation and 
mesons in RPA, we determined the two coupling constants separately by fitting the meson masses in the vacuum, and 
then investigated the phase structure for chiral symmetry restoration and pion superfluidity, especially the isospin 
effect on the chiral phase transition line and the minimum isospin chemical potential for pion superfluidity 

The main conclusions are: 

1) The two chiral phase transition lines in the T — plane predicated by the NJL model without 17a (1) breaking 
term is cancelled by the strong Ua (1) breaking term at low isospin chemical potential in real world (N c = 3). Only 
when the isospin chemical potential \ij is close to the critical value \i\ of pion condensation, there is probably the 
two-line chiral structure, depending on the ^/-meson mass. Therefore, in relativistic heavy ion collisions where the 
typical iii value is much less than it looks impossible to realize the two-line structure. However, it can be realized 
in large N c QCD. 

2) The critical isospin chemical potential for pion condensation in NJL model is exactly the pion mass in the vacuum, 
,Uj = to^ independent of the model parameters, the regularization scheme, the ?7a(1) breaking term and the color 
degree of freedom N c . 

The temperature and baryon and isospin density dependence of the chiral and pion condensates as well as the 
meson masses, and the extension to flavor SU3 NJL model are under way. 
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